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ABSTRACT 

We upgraded the chemical network from the UMIST Database for Astrochemistry 2006 to include isotopes such as 13 C and ls O. 
This includes all corresponding isotopologues, their chemical reactions and the properly scaled reaction rate coefficients. We study 
the fractionation behavior of astrochemically relevant species over a wide range of model parameters, relevant for modelling of 
photo-dissociation regions (PDRs). We separately analyze the fractionation of the local abundances, fractionation of the total column 
densities, and fractionation visible in the emission line ratios. We find that strong C + fractionation is possible in cool C + gas. Optical 
thickness as well as excitation effects produce intensity ratios between 40 and 400. The fractionation of CO in PDRs is significantly 
different from the diffuse interstellar medium. PDR model results never show a fractionation ratio of the CO column density larger 
than the elemental ratio. Isotope-selective photo-dissociation is always dominated by the isotope-selective chemistry in dense PDR 
gas. The fractionation of C, CH, CH + and HCO + is studied in detail, showing that the fractionation of C, CH and CH + is dominated by 
the fractionation of their parental species. The light hydrides chemically derive from C + , and, consequently, their fractionation state is 
coupled to that of C + . The fractionation of C is a mixed case depending on whether formation from CO or HCO + dominates. Ratios 
of the emission lines of [C II], [C I], 13 CO, and H 13 CO + provide individual diagnostics to the fractionation status of C + , C, and CO. 

Key words. Astrochemistry - ISM: abundances - ISM: structure - photon-dominated region (PDR) - ISM: clouds 



1. Introduction 

Astronomical observations of molecules and their respective iso- 
topologues reveal, that abundance ratios of the main species to 
their respective isotopologues may differ significantly from e.g. 
solar system isotope ratios. While isotopic fractionation in the 
interstellar medium is widely discussed in the framework of deu- 
terium chemistry, its relevance for the 12 C/ 13 C ratio in various 
species usually gains much less attention. In this paper we inves- 
tigate chemical fractionation in the context of models of photo- 
dissociation regions (PDR) with the focus on effects that result 
from introducing 13 C into the applied chemical network. 
The most important fractionation reaction is 



13 



C+ + CO 



C+ + 13 CO + AE = 35K 



(CI) 



(see I Woods & Willacvl [2009. and references therein). At high 
temperature back and forth reaction are equally probable and no 
apparent deviation from the elemental isotope ratio takes place. 
The lower the temperature gets, the less probable the back reac- 
tion becomes, resulting in a one-way channel shifting 13 C into 
13 CO and decreasing the abundance ratio of 12 CO/ 13 CO. At the 
same time the abund ance ratio of 12 C + / 13 C + is shifted oppositely. 
Lang er et al.l (Q~984) performed numerical calculations of a time- 
dependant chemical network for a variety of physical parameter 
and given kinetic temperatures concluding that chemical frac- 
tionation of carbon bearing species is of increasing significance 
the lower the temperate is, confirming the relevance of zero- 
point energy differences of a few ten K at low temperatures. In 
dark cloud models, where radiation is usually neglected, the ki- 
netic temperature is the major parameter in opening and clos- 
ing reaction channels. For a given density, the chemical network 
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only depends on the temperature and the cosmic ray ionization 
rate (cr G an d history in case of time-dependant calculations). 

Carbon fractionation in mole cular and diffuse c loud s has 
been discussed sy stematically by Keene et af] ([1998) and iLisztl 
(2007). Following Wakelam & Herbst (2008) we assume a stan- 
dard elemental abundance of 13 C/ 12 C of 67 in the solar neigh- 
bourhood. Observi ng 13 [Cl] and 13 C ls O in the Orion Bar 
iKeene et al.l d!998l) found little evidence for chemical frac- 
tionation. Their observations showed a slight enhancement of 
c 18q/13 C 18q =75 and no enhancement of [ 13 CI]/[CI] rela- 
tive to the standard elemental abundance while the chemical 
models predicted the opposite. The systematic study of the 
c i8 /13 C isq bv lLanger & Penziasl (fT990l[T99l showed a sys- 
tematic gradient with Galactocentric radius, i.e. significantly 
higher 13 C abundances in the inner Milky Way, but also varia- 
tions b etween 57 and 78 at the solar circle. Wou terloot & Brandl 
(119961) showed that the trend continues to the outer Galaxy 
with ratios above 100 in WB89-437. Optical spectroscopy of 
13 CH+ in diffuse clouds has shown that the 13 CH+/CH+ ra- 
tio matches the elemental abundance ratio in t h e solar vicin- 
ity very closely (see e.g. ICenturion et al.1 0*9951) . ILisztl (120071) 
showed that the fractionation reaction ( IC U is even the domi- 
nating CO destruction process for high and moderate densities. 
They also report fractionation ratios of the CO column density 
between 15 and 170 with a tendency for the ratio to drop with 
increasing column density. Their analysis is based on observa- 
tion of diffuse clouds and covers a CO column density range 
up to iV( 12 CO) » 2 x 10 16 cm- 2 and densities < 100 cm" 3 . 
Here we complement this study by concentrating on PDRs with 
higher densities (n > 10 3 cm -3 ). 

This paper is organized as follows: in Sect.|2]we will briefly 
overview the KOSMA-r PDR model which has been used to per- 
form the model computations. The updated isotope chemistry is 
described in detail in Sect. 12.21 In Sect. [3] we will present results 
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from our model calculations. In a second paper (Ossenkop f" et al.1 
I2012L Paper II) we present observations of the [Cll]/[ 13 C II] ra- 
tio in various PDRs and investigate the fractionation ratio of C + 
in more detail. 



2. The KOSMA-r PDR model 

A large number of numerical PDR codes is presently in use 
and an overview of many established PDR models is pre- 
sented inlRollig et all d2007lf T l We use t he KOSMA-r PDR code 
dStorzer et alj 11996; Rollig et al. 20060 to numerically solve 
the coupled equations of energy balance (heating and cooling), 
chemical equilibrium, and radiative transfer. The main features 
of the KOSMA-t PDR model are: a) spherical model symme- 
try, i.e finite model clouds, b) modular chemistry, which means, 
that chemical species can easily be added or removed from 
the network and the network will rebuild dynamically, c) iso- 
tope chemistry including 13 C and ls O, and d) optimization to- 
ward large model grids in parameter space, allowing for ex- 
ample to build up any composition of individual clouds in or- 
der to simulate clumpy material (for details see [Cubick et alJ 
2008). The KOSMA-t results can be accessed on-line at: 
|http : //www . astro . uni-koeln . de/~pdr| 

2.1. Model physics 

Individual PDR-clumps are characterized by the total gas den- 
sity n at the cloud surface, the clump mass M in units of the 
solar mass, the incident, isotropic far ultraviolet (FUV: 6 eV < 
E < 13.6 eV) intensity x, given in units of the mean interstel- 
lar radiation field of Draine (1978), and the metallicity Z. We 
assume a density power-law profile n(r) = no(r/R to t)~ a for 
-Rcore < r < i£tot> an d n(r) = const, for < r < i? CO rc- The 
standard parameters are: a — 1.5, i? C orc = 0.2i? to t, roughly 
approximating the structure of Bonnor-Ebert spheres. Figure Q] 
shows the applied density structure. In Appendix ID1 we describe 
how to compute mean column densities for spherical model re- 
sults. 

Excitation of the H 2 molecule is computed by collapsing 
all rotational levels with the same vibrational quantum number 
into a corresponding, virtual v level. We then solve the detailed 
population problem accounting for 15 ground-state levels {v = 
— 14) and 24 levels from the Lyman band as well as 10 Werner 
band levels. We assume that chemical reactions with the popula- 
tion of vibrationally excited H 2 have no activation energy barrier 
to overcome. This is especiall y important for specie s such as CH 
(iRollig et alJl2007l) and CH + (lAgundez et alJl2010h . The heating 
by collisional de-excitation of vibrationally excited H 2 is cal- 
culated from the detailed lev el population. Photo-ele ctric heat- 
ing is calculated according to lBakes & Tielensl (fl994). Overall, 
we account for 20 heating and coo ling processes. For a de tailed 
description see lRollig et ail f2006) and iRollig etal] d2007l) . All 
model results in this paper are for single-clump models without 
contribution from interclump gas and without clump superposi- 
tion. 



The comparison results and code descriptions are available at 
|http : //www .astro . uni-koeln . de/pdr- comparison! 

2 We recently updated KOSMA-r to self-consistently account 
for various dust compositions and dust size distributions including 
wavelength-dependent continuum radiative transfer, dust temperature 
compu tation and photo-electric heating. For details see Rolli g et al] 
d2012[) . However, the update was not fully available when we started 
this study, so we do not fully use the updated code capabilities here. 




Fig. 1. Radial density profile of a model clump. The dashed lines 
show the density profile along various lines of sight (for different 
impact parameters p). 



2.2. Model chemistry 

In KOSMA-t we solve pure gas-phase steady-state 
chemistry with the except ion of H 2 forming on grains 
([Sternberg & Dalgarnol I"l995). It relies on the availability of 
comprehensive databases of chemical reaction rate coefficients. 
Today a few databas es are publicly available, such as UDfAEI 
(IWoodall et al.ll2007|. the Ohio database OSlfl and the NIST 
Chemistry Webboolfl There are also efforts to pool all available 
reaction data into a unified database (KIDA: Kinetic Database 
for Astrochemistryfl 

In the following we use UDfA06. Retrieved at the 
09/29/2009 the database consists of 4556 reaction rates, involv- 
ing 421 chemical species (420 species + electrons). 34 of these 
reactions are pr esent with 2 o r 3 entries, valid in different tem- 
perature ranges (R611igl l2011h . All species in UDfA06 are com- 
posed of the main isotopes only. 

2.2.1. Isotopization 

In order to calculate chemical reactions involving different iso- 
topologues, i.e. molecules that differ only by their isotopic com- 
position, it is necessary to extend the chemical reaction set. For 
example, reaction 1 in UDfA06: 

H + CH — > C + H 2 

becomes 

H + 12 CH — > 12 C + H 2 
H + 13 CH — > 13 C + H 2 

We developed a software routine to automatically implement 
isotopic reactions into a given reaction sefl. A similar automatic 
proced ure was used by iLe Bourlot et alj d!993l) : iLe Petit et al.l 
(2006). We applied our routine to the UDfA06 reaction set, but 

3 |http : / / www . udf a ■ net| 

4 |http : / /www . physics ■ ohio- state . edu/~eric/researchTh" 

5 |http : / /webbook . nist ■ gov/chemistry/| 

6 |http : / /kida . obs . u-bordeauxl . f r"71 

7 We realized the isotopization routine in Mathematica©by Wolfram 
Research. 
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it can be applied to any given set of chemical reactions. The rou- 
tine features are: 

- inclusion of a single 13 C and a single I8 isotope (multiple 
isotopizations are neglected in this study) 

- UDf A often does not give structural information, for instance 
C 2 H 3 does not distinguish between linear and circular con- 
figurations (1— C 2 H 3 and c— C 2 H 3 ). In such cases we con- 
sider all carbon atoms (denoted by C n ) as indistinguishable. 
However, if structure information is provided we account for 
each possible isotopologue individually: 

- molecular symmetries are preserved, i.e. NC 13 CN = 
N 13 CCN, but HC I8 OOH ^ HC0 18 OH 

- functional groups like CH n are preserved (see also 
IWoods & Willacvll2009h 

- when the above assumptions are in conflict to each other we 
assume minimal scrambling, i.e. we choose reactions such, 
that the fewest possible number of particles switch partners. 

- we favor proton/H transfer over transfer of heavier atoms 

- we favor destruction of weaker bonds 

In Appendix [A] we describe in detail how isotopologues were 
introduced into the chemical network. The rescaling of the newly 
introduced reaction rates is described in Appendix lB^I 

2.2.2. Choice of the chemical data set 

Some data sets, e.g. OSU, are more relevant for the cold ISM 
while others, for instance by including reactions with higher ac- 
tivation energy barriers, are better suited to describe the warm 
ISM. Consequently, the results of chemical model calculations 
may differ significantly depending on the applied chemical data 
set. Reactions that can be found in different chemical sets can 
have significantly different rate coefficients among the vari- 
ous sets. Even very prominent reactions, such as the photo- 
dissociation of CO (from now on we will omit the isotopic su- 
perscript when denoting the main isotope) are listed with very 
different rate coefficients: UDfA06 gives an unshielded rate co- 
efficient of 2 x 1CT 10 s~\ OSU and KIDA give 3.1 x 1CT 11 s _1 . 
This is a huge difference and will lead to a significantly different 
chemical structure of a model PDR. It is not the purpose of this 
paper to perform a detailed analysis of how the choice of a chem- 
ical data set affects PDR model results. However, we show in 
AppendixIClhow the (isotope-free) chemistry of the main species 
discussed later in this paper changes for three different che mical 
sets. For a similar discussion see also Wakel am et al.l (120121) . 

2.2.3. Isotope exchange reactions 

In this frame reaction (IC II ) turns into two reactions, one for the 
forward, one for the back reaction, where both have the value of 
a (the rate coefficient for reaction (IC U is fc j^jy^ = 4.42 x 
10 -iO(7y3ooiq- - 2 9), but where the back rea ction is sup- 
presse d by the factor exp(— 7/T) with 7 = 35 K. IWatson et al.l 
(1 1 976 ) proposed that carbon isotope transfer between interstellar 
species can occur as a result of reaction ( IC II ). At low tempera- 
tures, reaction ( IC Q transfers 13 C isotopes from 13 C + — > 13 CO, 
enhancing the abundances of 13 CO and 12 C + . This reaction needs 
three main ingredients: sufficient amounts of 13 C + and CO and 
temperatures well below 100 K. In the PDR context, reaction 

The isotopized chemical data set as well 
as the isotopization routine are available online: 

|http : //www .astro . uni-koeln . de/kosma-tau| 



(IC U is especially interesting, because in the outer transition 
regions, where CO is still strongly dissociated, the 13 C + abun- 
dance is very large while I3 CO is very rare. Even small num- 
bers of I3 CO products from reaction ( IC U will have a signif- 
icant influence on the total 13 CO abundance and thus on the 
[CO]/[ 13 CO] ratio. Vice versa, 13 C + will be depleted strongly, in- 
creasing rC + l/r 13 C + l. [Smith & Adaml ll 19801) measured another 
isotope-exchange reaction: 

HCO+ + 13 CO H 13 CO+ + CO + AE = 9 K (C2) 

with less effect on the hotter parts of the PDR, due to the low 
differe nces in back and for ward reaction rates at higher temper- 
atures. lLanger et al.l (1 19841) tabulated reaction rates and reaction 
enthalpies for the various isotopic variants of reactions dC 11 ) and 
(IC 2\ and we use their values in our calculation s. Slightly d iffer- 
ent reaction rate coefficients are also given bv iLiszt) (120071) and 
IWoods & Willacy! ([2009) but the differences are small. 

3. Application 

3.1. Model parameter grid 

We test the outcome of the isotopic network under various condi- 
tions by computing a large grid of models spanning the possible 
parameter space. Our chemical network consists of 198 species, 
involved in a total of 3250 reactions. We did not include 18 into 
the chemistry here. 

We separate the model parameters into two sets: fixed and 
variable. The fixed parameters determine the fundamental phys- 
ical and chemical conditions for the model clouds, e.g. gas den- 
sity profile parameters a and i? CO rc/ Rtot, cosmic ray ioniza- 
tion rate (cr of molecular hydrogen elemental abundances Xi, 
metallicity, and dust composition. The variable parameters com- 
pose the final model parameter grid. A common set of vari- 
able parameters is: total surface gas density no = n(R to t) = 
fin + 2nn 2 , cloud mass M, and ambient FUV field strength 
X in units of the Draine field dDraindll978l) . For a given den- 
sity law, a and f c — R CO i C / Rtot, the total cloud radius i? to t = 
5.3 x 10 18 y/M/n cm and the maximum (radial) column density 
-Wmax = 4.7rai?cm -2 (see also App.lDl. 

For the present study we varied the clump parameters no, 
M, x and kept all other parameters constant. Table Q] gives an 
overview over th e used parameters. We assume a dust composi- 
tion according to Weingartne r & Draind d2001l) (entry 7 in their 
Tab. 1, which is equivalent to i?v = A.\ /E-q-v = 3.1). From 
the extinction cross section of each dust component we com- 
pute an average, effective FUV dust cross section per H <td- For 
a given total gas column density AT tot , then follows: rpuv = 
Ntot&T) > an d Ay = <7D-/Vtotl-086/3.08. The term in the de- 
nominator corrects from visual to FUV extinction. Note, that the 
elemental abundances of carbon show an elemen tal ratio (ER) 
of 67 , close to the average ratio in the local ISM (Shefferet al. 
120071) . We computed 168 models. The computation times per 
model range from 36 to 930 minutes with a median of 100 min- 
utes 

3.2. Structure of the C+/C/CO transition 

The chemical stratification of C + -C-CO, the carbon tran- 
sition (CT), is a well known signature of PDR chemistry 
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Table 1. Overview of the most important model parameter. All 
abundances are given with respect to the total H abundance. The 
numbers in parentheses indicate powers of ten. 



Model Parameters 



He/H 


0.0851 


Asplund et al. (2005) 


O/H 


2.56(-4) 


Wakelam & Herbst (2008) 


C/H 


1.2(-4) 


Wakelam & Herbst (2008) 


13 C/H 


1.8(-6) 




N/H 


6.03(-5) 


Asplund et al. (2005) 


S/H 


3.5(-6) 


Goicoechea et al. (2006) 


Z 


1 


solar metallicity 


CCR 


5(-17) s _1 


CR ionization rate 


Rv 


3.1 


Weingartner & Draine (2001) 




1.75(-21) cm 2 


UV dust cross section per H 


< A(\)/A v > 


3.339 


mean FUV extinction 


TUV 


3.074^ v 


FUV dust attenuation 


Vb 


1 kms" 1 


Doppler width 


no 


10 3,...,6 cm -3 


total surface gas density 


M 


10 -2,...,3 Mq 


cloud mass 


X 


10 0...,6 


FUV intensity w.r.t. 






Draine (1978) field^, 


a 


1.5 


density power law index 




0.2fltot 


size of const, density core 


Ntat/Av 


1.62(21) cm -2 






Av 

Fig. 2. CO shielding factorsT^of a model clump with the follow- 
ing model parameters: uq — 10 5 cm~ 3 , M — 100 Mq, x = 10 
(CO: solid line, 13 CO: dashed line). 

(e.g. iHollenbach & Tielensl [l999). The stratification is the re- 
sult of a number of competing formation and destruction pro- 
cesses like photo-dissociation, dissociative recombination, and 
others. Photo-dissociation of CO by FUV photons is a line 
absorption process, and as such subject to shielding effects 
( Ivan Dishoeck & Blackl 1 19881) . Numerically, this can be de- 
scribed by shielding factors that depend on the column den- 
sities of dust and of all species that absorb at the frequency 
of the line. The shielding factor multiplicatively enters the 
photo-dissociation rate, i.e. describes the reduction of the photo- 
dissociation through the line absorption. In the case of carbon 



Because of the isotropic illumination one has to compute the aver- 
age shielding over 4n at each Av- Optical thickness and scattering then 
leads to a mean shielding factor at the cloud surface of approximately 
0.5. 
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Fig. 3. Chemical structure of a model clump with the following 
model parameters: n = 10 5 cm~ 3 , M = 100 M , x = 10. 
(main isotopologue: solid line, 13 C isotopologue multiplied by 
ER=67: dashed line). 



mono xide the shielding depends on the columns o f H 9 and CO 
(e.g. Ivan Dishoeck & Blackl Il988t iLee et al.lll996t IWarin et al.l 
1996: IVisser et al.l2009H 13 r The self-shielding of the CO leads to 
stronger photo-dissociation of the rarer isotopologues at a given 
AvO This is also shown in Fig. [2] Because of the lower column 
density of I3 CO, with respect to the main isotopologue, it takes 
a larger cloud depth for 13 CO to become optically thick than it 
takes for CO. From Fig.|2]it can be seen that photo-dissociation 
of 13 CO is still strong at Ay w 0.3 where CO is already optically 
thick. 

The physical conditions, such as density structure and FUV 
illumination, determine where the CT is situated. For the pur- 
pose of this paper we define CT as the position in a cloud where 
n(C + ) = n (CO). The details of the C + -C-CO structure are 
changed by various effects. C + remains the least affected species. 
At the outside of the cloud photo-ionization turns basically all 
carbon into C + . The strength of the FUV field and the attenu- 
ation by dust determine the depth where recombination domi- 
nates over ionization and the C + abundances decreases. Carbon 
now becomes distributed between numerous species, but once 
shielding of CO becomes effective, usually at Ay > 1, the large 
majority of all carbon atoms is bound into carbon monoxide. 
Both species, C + and CO, are quite insensitive to changes in the 
chemistry or the temperature structure, at least in cloud regions 
where they dominate the carbon population. However, in regions 
where they represent only a minor fraction of all carbon species, 
their chemical structure may depend sensitively on details of the 
cloud chemistry and physics. For example, increasing the H 2 for- 
mation efficiency on hot dust grains, i.e. at low Ay, increases the 
corresponding H 2 formation heating efficiency, leading to higher 
gas temperatures in these cloud parts. This can produce an in- 
crease of the CO population in the hot gas a nd consequently 
produce strong emissio n of high- J emission dLe Bourlot et al.l 
2012; Rol lig et al.ll2012l) . The same effect can also positively af- 
fect the abundance of light hydrates, such as CH + and CH, which 
are primarily formed in these regions. 

13 Th e most recent set of shielding rates is pr edicted by I Visser et all 
(2009). In this work we still used the rates from van Dishoeck & Black 
( 1988). Under most parameter conditions the model results differ only 
marginally (~ 10%) when switching to the new shielding rates. 

14 12 CO could also shield its less abundant relatives, if their absorp- 
tion lines are sufficiently close together. However, only very few of the 
efficiently dissociating lines overlap, too few for mutual shielding to be 
important (War in et aljl996h . 
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Atomic carbon is the species that is probably most affected 
by changes in the chemistry and physics, because it is involved 
in the chemistry of C + as well as CO and chemically consti- 
tutes a transitional species. It is the major carbon species, that is 
least understood. Observations and model predictions of the spa- 
tial distribution of C show big differences. The classical C + -C- 
CO stratification, with C being sandwiched between its two big 
brothers is hardly observed at all. Instead, atomic carbon shows 
a widespread distribution that remains to be understood (e.g . 
Kramer et alJl2004t iMookeriea et al.ll2006t iKramer et all feOQS: 
Rolli g et alj|201 UlMookeriea et alj|2012l) . 

3.3. Cloud fractionation structure 

Unfortunately it is impossible to discuss the fractionation for all 
species from our chemical dataset so that we focus on a few 
molecules of particular astronomical interest. A complete cover- 
age of the fractionation ratio (FR) for the selected species in our 
model grid is presented in Appendix[E] 

3.3.1. C + and CO 

Figure[3]shows the chemical structure of the main carbon species 
in the model clump. Solid lines show the main isotopologues, 
dashed lines the 13 C variants. In the outer parts of the clump, 
most of the carbon is in the form of C + . The CT for the model 
shown in Fig.[3]is at Ay w 0.2 where the CO photo-dissociation 
rate has dropped sufficiently in order to build large quantities of 
CO, which leads to a steep decline in n(C + ) and n(C). In Fig. [4] 
we show the corresponding FR of the species from Fig. [3] We 




Fig. 4. Fractionation structure of the same model clump as 
shown in Fig. [3](solid lines: fractionation ratios n(X)/n( 13 X), 
dashed line: kinetic gas temperature (right axis). 

note a number of features: 

1. the FR of C + is always larger than or equal to the ER., i.e. 
13 C + is always under-abundant with respect to C + 

2. the FR of C + equals the ER at low Ay 

3. the FR of C + increases significantly at large Ay 

4. the FR of CO is always smaller or equal to the ER except for 
conditions described in Sect. 13.3.21 

5. the FR of CO deviates the strongest from the ER at low Ay 
and equals the ER at large Ay 

6. the FR of C and HCO + show mixed behavior 

Points 1-3 are a direct consequence of reaction ( 1C 11 1. Any 
fractionation of C + has to be a direct result of this reaction. C + 



stands at the beginning of the chemical chains and there is no 
other direct mechanism acting in the opposite direction, such as 
isotope-selective photo-destruction. At low Ay temperatures are 
sufficiently large for the reaction to work equally well in both 
directions. Once the temperature drops below 50 K the back re- 
action becomes less probable and the FR increases rapidly and 
can be kept at large values as long as enough 13 C + ions are avail- 
able to feed the reaction. Deep inside the clump, the FR is a 
result of the balance between 13 C + formation via He + + 13 CO 
and destruction via the fractionation reaction. Dominance of re- 
action ( IC Q automatically leads to fractionation, in this case to 
C + enrichment relative to 13 C + . The absolute magnitude of the 
FR is controlled by the He + abundance which is a direct result 
of the cosmic ray ionization rate. In Fig. [4] this can be seen by 
the roughly constant FR of C + deep inside the cloud. The same 
qualitative behavior is visible for all other model parameters in 
our model grid. Only the particular position, width and height of 
the FR peak of C + varie s with density and FUV field strength. 
IWoods & Willacvl (|2009) find the same fractionation behavior in 
their protoplanetary disk model calculations. At the surface C + 
is not fractionated while at large depths they find FR>ER. 

The fractionation of CO is different because a second 
isotope-selective process is at work, the shielding of CO and 
13 CO from FUV photons. We would expect that if photo- 
dissociation was the dominant process, i.e. if 13 CO photo- 
dissociation was relatively stronger than the photo-dissociation 
of the main isotopologue, then this would result in FR(CO) > 
ER. This is not the case for any model clump in our calculations. 
It happens in th inner clouds with n < 10 2 cm~ 3 as discussed 
by iLisztl (120071) . All models show a FR{CO) < ER indicating 
that the FR is dominated by the chemistry, i.e. by reaction ( IC II ) 
which can only produce FR < ER. Exceptions to this behavior 
are discussed in Sect. 13.3.21 




Av 



Fig. 5. Reaction rates for the same model clump as shown in 
Figure [3] (photo-dissociation of CO: black, solid; 13 CO: black, 
dashed; back and forth reaction rates in reaction dC II ): red). 



This is also visible from Fig. [5] where the destruction of 
CO (black, solid line) and 13 CO (black, dashed line) via photo- 
dissociation is compared to the respective formation via reac- 
tion ( IC II ) (red lines). For CO photo-dissociation is the major 
destruction process until Ay > 0.1. Then chemical destruc- 
tion by the fractionation reaction takes over. Formation via the 
fractionation reaction is not the dominant formation channel of 
CO for most of the clump. Only for the small Ay range of 
0.1 < Ay < 0.6 CO formation is dominated by the fractionation 
reaction. This is different for 13 CO where destruction via photo- 
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dissociation is weaker than formation by reaction (IC U through- 
out the clump. Both, formation and destruction of 13 CO is gov- 
erned by reaction ( IC II ) (both red lines in Fig. [5J. At Ay > 2 
electron recombination with H 13 CO + becomes the main forma- 
tion channel. Hence, for the whole low Ay part of the clump, 
the I3 CO abundance is controlled by the chemical fractionation 
reaction and accordingly, I3 CO is significantly enriched relative 
to CO. 

Figure [6] compares the FR of C + (left panel) and CO (right 
panel) for n = 10 3 cm -3 and M — 1 M Q . Each sector in the 
figure corresponds to a different FUV intensity \- The FR, as 
function of the relative clump radius r/R tot , is color coded, 
ratios within ±10% of the ER are shown in green. Figure [6] 
gives a visual summary of the analysis above. The pie chart 
representation visualizes the relative contribution of the abun- 
dance profile at different radii of our spherical clumps to the 
integrated clump ratio. At n = 10 3 cm~ 3 , fractionation of C + 
only occurs deep inside the clump, while CO shows most frac- 
tionation further out. However, we also note, that for higher 
values of x, the cloud is so hot that no fractionation occurs 
any more. The stronger the FUV intensity the deeper the domi- 
nance of the photo-dissociation of CO. No shielding or selective 
photo-dissociation can yet take place and the FR equals the ER. 
However, in case of n — 10 3 cm" 3 and \ > 10 2 > CO can not 
be shielded efficiently and most of the carbon is locked in its 
ionized form. This is different for higher densities. 

The situation at higher gas density is presented in Fig. [7] 
The fractionation of C + is much more prominent and dominates 
a much larger clump volume compared to lower densities. As 
explained above, cold (T<100 K) C + is always fractionated with 
FR^ER. This is true for the entire parameter grid. CO on the 
other hand, requires CO and sufficient 13 C + to become fraction- 
ated. These conditions are only met in a limited radius range, 
that is pushed to larger depths if \ increases. Deeper inside and 
further outside, the FR of CO equals the ER. 

The effect of different clump masses is easier to understand. 
Adding mass is effectively equivalent to adding shielded, cold 
material to the clump, since it approximately requires a constant 
column of gas to attenuate the FUV radiation. Once this column 
is reached any additional material will be shielded and therefore 
located in the center of the clump. The appendix gives sector 
plots such as Figs. [6] and [7] for the full parameter grid. 

As an additional effect, the 13 CO recombination occurs at 
lower values of Ay than the recombination of the main isotopo- 
logues, despite the lower shielding capabilities of 13 CO com- 
pared to CO. This is true for all models in our parameter grid. 
Across our parameter grid, the CT of the main isotopologue oc- 
curs at a log(iVco) = 15.8 ± 0.4, while for the 13 C variant 
( 13 CT) log(AT co ) = 15.4 ± 0.4. The difference between both 
is smallest for models where photo-dissociation is more impor- 
tant, i.e. for larger values of x- The comparison of the CO col- 
umn density at the CTs of CO and 13 CO for all models in our 
parameter grid is shown in Fig. [8] The colors of the data points 
represent the FUV field strength x of the respective model. 

The resulting impact on the column densities is shown 
in Fig. [9] Each symbol represents the column density 
(iV(CO))/(iV( 13 CO)} of a model clump with given density n, 
mass M, and FUV irradiation x (clump column densities are 
defined in Appendix ID1 Clumps with low CO column densities 
deviate most from the ER. Measurements of the column den- 
sity ratio N(CO)/N( 13 CO) are usually performed on diffuse or 
translucent clouds and thus naturally confined to a low TV (CO) 
regime. 




0.05 0.10 0.50 1.00 5.00 

Ay 12 CT 



Fig. 8. Comparison of the Ay i2 CT and Ay i3ct> i- e - the Ay t 
where u(C+) = n(CO) and n( 13 C+) = n( i3 CO) for all mod- 
els in our parameter grid. Each data point is colored accord- 
ing to the ambient FUV field x- The solid line corresponds to 

A v ^i2 CT — A v ,i3ct- 

The low column density region of Fig.|9] is roughly consistent 
with U V absorption-line observations by ISonnentrucker et all 
(120071) . For translucent clouds they found an anti-correlation 
of (7V(CO))/(N( 13 CO)) with (N(CO) \ in t he ra nge of 
10 14 cm" 2 < (N(CO)) < 10 16 5 cm" 2 . iLisztl (120071) found 
15 < N(CO)/N( 13 CO) < 170 with a tendency for the ratio 
to decline for higher column densities and a total CO column 
densities of a few 10 16 cm -2 from Galactic CO absorption and 
emission at 1.3 and 2.1 mm wavelengths f or clouds with a tota l 
CO column density AT(CO) < 10 16 cm~ 2 . IShefferetal] (f2007) 
showed, that UV data toward diffuse/translucent lines of sight 
can give 0.5 < FR{CO) /ER < 2. 

Clumps with stronger FUV fields show almost no fractiona- 
tion, either because the molecular inner parts are so small that the 
gas temperatures throughout the clump are too high or because 
fractionation only affects the CO at the outer clump regions but 
not the bulk of the CO gas. For sufficiently large CO column 
densities, the column density ratio of CO/ 13 CO turns out to be a 
relatively good tracer of the elemental abundance ratio of a given 
cloud. 

The fractionation of the column density ratio of C + is shown 
in Fig. [10] The lower FUV models show the largest FR while the 
models with very large FUV fields have a FR=ER. Models with 
low density show a much weaker fractionation (see also Fig. IE. Il l 
because the molecular part of the low density clouds contributes 
less to the total C + column density. Increasing the model den- 
sity (larger symbols in the figure) for given model mass and 
FUV illumination moves the models in Fig. QJJto the top left 
because larger densities lead to a stronger shielding of the gas 
from the FUV and therefore a CT closer to the clump surface. 
Consequently, the total C + column density is decreased. On the 
other hand, the larger amount of cold CO gas acts in favor of 
the stronger C + fractionation. However, these strongly fraction- 
ated model clumps show the lowest C + column density making 
observations difficult. Observability will be further discussed in 
Sect. 133] 

As an additional consequence, when keeping density and 
FUV constant, the FR of C + is proportional to the clump mass. 
Increasing the clump mass leads to an increased FR, moving 
model points, visible as different symbols in Fig. [10] to the top 
right. 
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Fig. 6. Fractionation structure as function of relative clump radius r/R tot for n = 10 3 cm -3 and M = 1 M©. Each sector 
corresponds to a different \ value. The FR is color coded, ratios within ±10% of the ER are shown in green. Left panel: FR of C + , 
the color scale goes from to 2000. Right panel: FR of CO, the color scale goes from to 100. 
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Fig. 7. Same as Fig.|6]for n = 10 5 cm 3 and M = 1 M . 



3.3.2. High density conditions with FR{CO)>ER 

For typical molecular cloud conditions isotope-selective photo- 
dissociation is never strong enough to keep 13 CO photo- 
dissociated while CO is already recombined. If isotope-selective 
photo-dissociation was dominating the CO chemistry one would 
obtain FR > ER, like it is found e.g. in diffuse clouds iLisztl 
(2007). For special conditions, particularly for densities of n > 
10 6 and sufficient FUV illumination, FR>ER is still possible 
across a limited Ay range. Figures [TT1 and [121 show the density 
and fractionation structure of such a clump. CO shows a first, 
relatively strong abundance peak at Ay — 0.6, before the CT. 
The gas temperature at this Ay is around 600 K. At this part of 
the clump, the FR of CO is higher than the ER. Closer examina- 
tion of the parameter reveals that similar behavior can be found 
for models with n > 10 6 and \ > 10 4 (compare Fig. IE. 3b . 

Detailed chemical analysis of these models shows that 
FR>ER occurs at significantly lower cloud depth than the CT 
(see Fig. [TTJ. The already mentioned CO abundance peak re- 
sults from dissociative recombination of HCO + , which is effec- 



tive in hot gas (ISternberg & Dalgarnolll989l) and gives rise to 
the CO peak before the CT At these cloud depths HCO + it- 
self is primarily formed by collisions of HOC 1 " and CO + with 
H 2 (CO + forms through C + + OH 
two main formation reactions: C + 



CO + 



HOC + + H). 



CO + 
H 2 



H and HOC + has 
> HOC + + H and 



The maximum CO abundance in this peak strongly increases 
with the total gas density. Even so, we have not yet reached the 
CT, the total CO column density from the cloud's edge to that 
peak position can already re ach values > 10 15 cm" 2 , wh ich is 
sufficient for self-shielding (Ivan Dishoeck & Blacklll988). Any 
CO self-shielding will be stronger for the main isotopologue 
than for 13 CO and if photo-dissociation is the main destruction 
process for both isotopologues than the stronger shielding of CO 
can give rise to FR>ER. If the destruction of 13 CO is controlled 
by reaction ( ICTT ) than FR<ER. 

To understand which conditions can lead to FR>ER we bal- 
ance the main formation and destruction processes of 13 CO. 
We already noted that formation via dissociative recombina- 
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Fig. 9. Scatter plot of the mean column density fractionation ratio vs. the mean column density of the main isotopologue of CO/ 13 CO 
of the whole parameter space. The model parameters n, M, and x> are coded as size, shape, and color of the respective symbols. 
The red line denotes the model elemental abundance [ 12 C]/[ 13 C]=67. 
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Fig. 10. Same as Fig.|9]for C + . 



tion of HCO + is a general requirement, otherwise CO self- 
shielding will not be effective. A second possible formation 
route is by the isotope exchange reaction ( IC II ). Destruction can 
be either via the back reaction ( IC II ) or by photo-dissociation. 
The gas temperatures are high enough to neglect the energy 
barrier of the back reaction and assume fc j^jTJ ^i s» k ^Q 

The rate coefficient for dissociative recombination of HCO + is 



fe(DR) = 2.4 x lQ- 7 {T/2,Q0K)-°- 29 s" 1 cm 3 , and C13 is the 
photo-dissociation rate of 13 CO. Hence 



fcfCT^,«(C + )n( 13 CO) + ? i( 13 CO)Ci3 = 
fc ICTV-)"( 13c+ )"( CQ ) + «(H 13 CO+)n( e -)A; (DR) . (1) 
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Fig. 11. Fractionation structure of the very dense model clump. 
The model parameters are n — 10 7 cm -3 , M — 0.01 Mq, % — 
10 5 . (main isotopologue: solid lines, 13 C isotopologue multi- 
plied by ER=67: dashed lines) 
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Fig. 12. Fractionation structure of the same model clump as 
shown in Figure QT| 



With ra(CO) = FR x n( 13 CO), n(e~) « n(C+) « 1.2 
10- 4 n tot , and n(C+) =£fix n( 13 C+) follows 



\ T u - 4 n("CO) n to t / 

(2) 

The last two terms in parentheses compete in changing the 
FR relative to ER. In the relevant cloud regime, each term lies 
between « 0.1 — 10, depending on the detailed conditions, and 
FR(CO) can increase to w 80 - 100. 

We emphasize, that this behavior is not just a sim ple compe- 
tition between chemistry and photo-dissociation as lLiszll (120071) 
described for diffuse clouds. A FR>ER in PDRs is the result of a 
local, chemically induced, dominance of the photo-dissociation 
of 13 CO over its chemical destruction. Only for 10 6 < n < 
10 7 cm~ 3 and 10 4 < n\ < 10 6 are the conditions such that in 
a narrow Ay range isotope-selective photo-dissociation leads to 
FR >ER. 



x=io< 



x=io 




Fig. 13. Same as Fig.|6]for C for n — 10 5 cm 3 and M = 1 M c . 



5 — 10 atomic carbon is dominantly formed through one of the 
following reactions: 



C+ +e" 

C+ + s 



c 

C + S+ 

c + o 



while destruction occurs mainly via photo-ionization. As dis- 
cussed, the FR of CO and C + behave oppositely in distinct cloud 
depths; FR(CO)<ER, mainly in outer layers, while FR(C + )>ER, 
somewhat deeper in. 

At low Ay both isotopic variants of atomic carbon will 
be dominantly formed from the recombination of their ionized 
forms and C will share the FR of C + . This is visible as first peak 
in FR(C).The weaker shielding of 13 CO and the strong fraction- 
ation of C + makes photo-dissociation of 13 CO the major forma- 
tion reaction for 13 C for Ay > 0.2. The isotope-selective shield- 
ing drives the FR of C toward lower values and gives rise to the 
dip in the curve in Fig. Qj] The magnitude of this diminishment 
depends on the differences in the shielding of 12 CO and 13 CO 
and the cloud depth where the shielding is still weak. If this dif- 
ference is still significant when FR(C) reaches its first peak it can 
push the FR to values smaller than the ER. 

For Ay > 0.2, until FUV shielding becomes strong, 13 CO + 
hv will be one order of magnitude faster than any of the other 
formation reactions while photo-dissociation of 12 CO will never 
be the dominant formation reaction. At Ay > 1, FR(CO)=ER 
and the FR(C) will increase again. Charge transfer between C + 
and S becomes the major formation reaction and consequently, 
C will share the FR of C + . At very large cloud depths cosmic- 
ray-induced photo-dissociation of CO becomes a main forma- 
tion reaction together with charge transfer between C + and SO 
for C and charge transfer between N + and 13 CO for 13 C. As a 
result, the FR(C) will slowly decrease with increasing Ay. 



3.3.3. C 

From Figs. |4] and [13] we see that the FR of C has a small regime 
at low Ay and low FUV intensities where FR<ER, while under 
all other conditions the FR>ER. The FR(C) starts to peak at the 
rise of FR(C + ) followed by a dip where CO turns to the ER and 
by a second peak at the declining flank of FR(C + ). Until Ay < 



3.3.4. CH + and CH 

Figure [14] summarizes the dominant formation and destruction 
channels of the discussed hydrocarbons. The arrows denote the 
primary reaction channels for CH + ,CH2, and CH3 (left panel) 
and for CH (right panel) across the model clumps. Figure [15] 
shows the FR of C + , CH + , CH£, CH^, and CH and Fig.QJshows 
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Fig. 14. Chemical network of the dominant formation and destruction channels for: left panel: CH + , CH 2 , and CHj; right panel: 
CH 



C+ ™£> C 2 N+ ^ C 2 ^ CS+ % HCS+ ^> CH 

From Fig.[15]and from the above chain of reactions it is obvious, 
that the fractionation ratio of C + will be handed down through 
the chain, unless other carbon species become involved. This is 
the case for CH. At very low Ay where the formation via C + 
H 2 is most important, and the FR is closely related to C which 
mostly equals the ER under these conditions. Once the FR of C + 
starts increase, it will affect the fractionation of all the related 
CH„ and consequently that of CH. 

Deeper in the cloud the same happens along a different 
chemical track. Recombination of CH3 and of HCS + (once 
Ay > 3) are the main formation reactions for CH. The chemical 
chain C + , C 2 N + , C 2 , CS + , HCS + , CH shares a common behavior 
of the FR (see also Fig. [T5J. As a side remark we would like 
to emphasize this chain as a good example of how the chemical 
networks of different elements, sulphur and nitrogen in this case, 
are mixed. Consequently it is important to include both networks 
to correctly compute the carbon chemistry. 

In Figs. Q~6] and [17] we show the FR of CH + and CH respec- 
tively. The FR(C + ) deviates from the ER only at larger values 
of Ay and column densities of species whose FR depends on 
FR(C + ) will only be affected if these deeper cloud regions con- 
tribute significantly to their total column density. For CH + this is 
the case for the same parameters where also C + is fractionated 
and Fig. [TBI shows their close chemical relationship. Deviations 
from the ER occur only for the low \ models with sufficient total 
column densities, similar to C + (see Fig. ITOt. 

The weak deviations of the FR(CH + ) from the ER is con- 
sistent with observations. ICenturion et all y995) found a mean 
value of the CH + / 13 CH + column density ratio o f 67 ±3 for five 
lines o f sight, very close to the interstellar ER. ICasassus et al.l 
(2005) report an average ratio of 78±2 from measurements 
along 9 lines of sig ht. Recent absorption-line observations by 
Ritch ev et alJ (1201 II) along 13 lines of sight through diffuse 
molecular clouds confirm a FR(CH + ) close to the ambient ER. 
They report total column densities of CH + of a few 10 13 cm -2 . 



the abundance profile of hydrocarbons and their respective iso- 
topologues. 

The reaction C + + H 2 + AE — > CH + + H requires an 
activation energy of AE=4600 K, but collisions with vibra- 
tionally excited H^ allow to overcome the barrier (Rolli g~et alJ 
120071: lAgundez et aTll2010l) . At the edge of the cloud, CH + is 
primarily produced from C + colliding with excited, molecu- 
lar hydrogen. At Ay ~ 10~ 3 , the proton exchange reaction 
C + + CH — > CH + + C together with ionization of CH become 
the main formation reactions. In those two regimes, the FR is 
controlled by C + (see also Fig. ll3Tl. At Ay ~ 1 the main forma- 
tion occurs via collisions of C with HCO + or H3 . 

Under very high density and very low x conditions, C + will 
be less abundant than C throughout the clump. As a result the 
dip between the two peaks visible in the FR{CH + ) in Fig. Q3] 
becomes much more prominent and can reach values below ER. 
This is visible in Fig. IE. 61 

The formation of CH originates at C + . Successive collisions 
with H 2 form the chain: 



CH 



cm 



CH 



At the end of the reaction chain dissociative recombination leads 
to CH and CH 2 . At low Ay, CH can also be formed from C, at 
high Ay it forms via dissociative recombination of HCS + . The 
reaction paths to CH (in order of cloud depths where they domi- 
nate) then are: 



C 



h.; 



C 



C+ ^>CH+' 




> CH -^U CH+ 




■> CH —} CH+ 
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Fig. 16. Same as Fig.|9]for CH + . 
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Fig. 17. Same as Fig.|9]for CH. 



CH remains abundant approximately until Ay approaches 
unity. As a consequence, a larger fraction of its total column 
density will be affected be the fractionation. This is visible in 
Fig- Ell Above a mean column density of 10 13 cm -2 , all mod- 
els show an enhanced FR. Because of the strong coupling to 
FR(C + ), weaker FUV models tend to have the strongest frac- 
tionation. Therefore, CH promises to be a good observational 
fractionation tracer as it combines enhanced FR with high col- 
umn densities. 



3.3.5. HCO + 

The fractionation of HCO + is special because it is affected 
by two processes acting in opposite directions. At very low 
A v HCO + is formed by H 2 collision with HOC + and CO + . 
Both precursors are not fractionated, thus FR(HCO + )~ ER (see 
Sect. 13.3.2b . A little deeper into the clump, the main formation 
reaction changes to CH + O — > HCO + + e~, thus its fractiona- 
tion indirectly depends on reaction ( IC 11 1. CH is strongly frac- 
tionated with FR > ER and passes down the fractionation to 
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HCO + /H 13 CO + 




HCO + . This fractionation peak is seen in Figs. l4l and [T9l With 
growing \ the peak is pushed to larger cloud depths. Once CO is 
sufficiently abundant, the reaction H3 + CO — > HCO + + H 2 
takes over as dominant formation reaction and the FR ap- 
proaches that of CO (compare with Fig. [7] right panel). At even 
higher values of Ay the gas temperature becomes very low 
(T < 10 K) and reaction ( IC 2\ starts to dominate formation and 
destruction of HCO + and pushes the FR significantly below the 
ER. In the appendix we show the FR of HCO + over a significant 
portion of our model grid. The central cloud regime with FR<ER 
is visible in all clumps that are sufficiently shielded from the ex- 
ternal FUV radiation. 
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Fig. 18. Same as Figure|3]for species HCO + , CH, CH + , CH 2 , and 



In Fig. [20] we show the column density fractionation ratio 
of HCO + . Fractionation of HCO + is strongest for clumps with 
large columns of cold gas, where reaction (IC 2b can contribute 
strongly to the total H 13 CO + abundance. Low mass and low den- 
sity models have a FR equal to the ER or slightly higher. For a 
given density, the FR is largely independent of the model mass, 
which is consistent with the HCO + and H 13 CO + density profiles 
shown in Fig. [18] which show an increase for large Ay and a 
roughly constant FR (see also Fig. 0J. Hence, the FR of HCO + 
is only marginally affected if the clump mass is increased. The 
model results show a strong correlation of the column density 
ratio with x- 



Fig. 19. Same as Fig. for HCO + for n = 10 5 cm" 3 and M = 
1M . 



3.4. Emission line ratios 

Even though column densities are no direct observables, they 
have to be derived from measured line strengths resulting 
from the full radiative transfer, including effects of a variable 
temperature of exciting collision partners and optical depths. 
However, the derivation of intensities requires numerous addi- 
tional assumptions, such as collision rates, details of the ge- 
ometry, assumptions on chemical and radiative pumping, and 
so on. Consequently, from a modellers perspective intensities 
have larger uncertainties than column densities. Therefore we 
only compute intensities for a selected subset of commonly ob- 
served species and transitions. Here, we directly compare the 
isoto pic ratio of clu mp averaged line intensities (for definition 
see lRollig et alj|2006l) which can be compared directly to obser- 
vations. We always assume unity beam filling factor. The large 
number of possible line combinations prohibits a complete pre- 
sentation. We give just a few examples to demonstrate that line 
ratios between various isotopologues can differ from the corre- 
sponding column density ratios. 

Fig. ED shows the ratio of T mb ([CII])/T mb ([ 13 CII]) with 

Tmb = / T m bdv (For 13 C + we summed over all hyperfine com- 
ponents.). While the column density ratio is larger than the ER 
for all model clumps the same is not true for the intensity ratio 
(IR). All models with \ > 100 show a IR<ER down to values 
of 38. Only models with very low x an d large densities have a 
IR>ER. For a given density and mass, the IR decreases with in- 
creasing x, because of the larger optical thickness of C + relative 
to 13 C + . [ 13 CII] observations of the Orion nebula bv lStacev et al.l 
<fl99lh show a very similar behavior and are consistent with our 
model predictions. They report line ratios between 36 and 122 
with the highest optical depths r = 3.31q 8 belonging to the 
lowest lin e ratio of 36 ± 9 and l ow optical depths where ratios 
are high. Bor eiko & Betzl d 19961) derive an intensity ratio of 46 
in M42, consistent with an intrinsic ER of 58^5 and an opti- 
cal depth of [CII] of 1.3. For a given x> the IR drops with de- 
creasing density and increasing mass. Increasing the mass for a 
given x an d n will add cool, molecular mass. Provided that C + 
is fractionated this will increase the FR as long as the [ 12 CII] 
line remains optically thin and decrease the FR once it becomes 
optically thick (compare Fig. [2TJ. 
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Fig.21. Same as Fig. for the intensity ratio T mb ([CII])/T mb ([ 13 CII]) vs.T mb ([CII]). 



In the previous section we showed, that the column density 
ratio of CO/ 13 CO is close to the ER of the clump for most of the 
parameter space. CO emission lines suffer much more from opti- 
cal thickness effects than most other species. This is also shown 
in Fig. [22] where we plot the ratio of the integrated intensities of 
COd-Oy^COd-O). For most of the models, this ratio is be- 
tween 1 and 10, much smaller than the ER. Quite a few models 
have intensities of the rarer isotopologue comparable to those of 
the main species. 



The low mass, high density models show an increasing IR 
with x> opposite to the high mass, low density models. This 
is a result of the CO abundance structure of these model (see 
Fig. fTTb . The lowest IR is shown at the highest model mass 
model (IR < 2) while the highest IR belongs to models with 
the lowest mass. Both can be explained as optical depth effect. 
Models with n > 10 7 cm~ 3 have an increasing T , m b(CO( 1 — 0)) 
with x and the IR increases for the lower mass models and de- 
creases for the higher mass models. The latter results from high 
optical thicknesses while the former results from the emission 
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Fig. 22. Same as Fig.EQfor T mb (CO(l 



0))/T mb ( 13 CO(l 



0))vs.T mb (CO(l-0)). 



of of very hot, strongly excited CO gas in the primarily ionized 
fraction of the clump. 

In Fig. [23] we plot the intensity ratio of CO lines 
r mb (CO(J ->• J - l))/T mb ( 13 CO(J ->• J - 1)) as a function 
of J for models with M — IOq and line intensities > 
. 1 K km s _ 1 . For the lowest transition, the ratio lies between 1 
and 10 and approaches ER for high values of J, when both lines 
become optically thin. The J where the ratio reaches the ER in- 
creases with n and \- For instance, models with n = 10 7 cm -3 
and x = 10 6 have a ratio of close to unity, until J > 10 and 
reaches ER at J > 25, while models with the same density and 
X = 10 3 reach ER already at J = 15. We find that either the CO 
lines are optically thick so that the IR is lowered below the ER or 
the 13 CO is too weak to be detectable. Only in dense clouds the 
ER is observable. The CO IR is therefore not a good diagnostics 
of carbon fractionation. 
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3.5. Diagnostics 

Finally, we studied several line ratios with respect to their di- 
agnostic value for the local FR of CO and C + as well as to the 
local ER. We already concluded that CH and HCO + appear to be 
sensitive tracers of the FR. Lacking collision rate coefficients for 
CH we only calculated HCO + (and H 13 CO + ) intensities for our 
model clouds. 

We selected transitions that are observable from the ground 
or through SOFIA (Stratospheric Observatory For Infrared 
Astronomy) and show at most moderate optical depths of 
a few. Figure [24] demonstrates how the emission line ratio 
Tmb([CII])/T mb (H 13 CO + (l - 0)) traces the column density 
FR of C + , C and CO. The figure reproduces the general trend 
of the FR as discussed in the previous sections. We find a clear 
dependence of the emission line ratios on the column density 
FR of C and C + , but not of CO. However, it is only partially 
applicable to observational data. All ratios above ~ 30000 are 
practically not observable as they correspond to H 13 CO + inten- 
sities below 10 mK km/s. If we exclude models with line inten- 
sities < 0.01 Kkms -1 fhe column density ratios of CO and C + 
show much weaker variations so that the emission line ratio only 
traces the fractionation of atomic carbon in the observable inten- 
sity range. We have repeated this exercise for many other line 
ratios that are available through ground based and satellite ob- 
servatories, and that are strong enough to be observabl43 All 
three panels in Fig. [24] show a transition around an emission 
ratio of 10 2 — 10 3 . This ratio corresponds to models with the 
highest T mb (H 13 CO+(l - 0)) « 1 - 2 Kkms^and FR(C) 
slightly below the ER. It requires x = 10 4 — 10 5 and densi- 
ties of n = 10 7 to reach these high values. Towards lower in- 
tensities, the models split into two branches: very high density 



Even though the PDR model takes optical thickness effects into 
account we tried, whenever possible, to find optically thin tracers in 
order to be as independent as possible from the local structure of the 
emitting source. 
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models with n = 10 7 cm -3 and decreasing FR(C), and mod- 
els with n = 10 6 cm -3 , \ > 100 and increasing FR(C). This 
parameter regime corresponds to model clumps with the high- 
est [CII] emission, increasing with n and x- Accordingly, the 
[CII]/H 13 CO + intensity ratio can increase to very high values if 
the density and the FUV field are high enough, albeit with the 
visible transition. 

Fig. [25] shows the result for the three ratios that seem to be 
suitable to trace the column density FR of C + , C, and CO individ- 
ually. The left panel shows T mb ( 13 CO(2 - l))/T mb ([CI] 610Aim ) 
on the x-axis and the CO column density ratio on the y-axis. The 
level energy of 13 CO (2-1) is 16 K, while the [CI] line requires 
24 K for excitation. Looking at the plot, we note, that any line 
ratio < 1 signals significantly fractionated CO column densities. 
Any emission line ratio higher than a few reflects a normal FR 
in the CO gas. All models with an emission line ratio smaller 
than w 1 have n < 10 4 cm~ 3 and will consequently host a small 
13 CO population but significant amounts of atomic carbon. 

The middle panel plots the column density fractionation ra- 
tio of atomic carbon versus T mb ([CII])/T mb (H 13 CO + (l - 0)) 
(same as in Fig. [24] but excluding all model lines too weak to 
be detectable). The model results show a strong correlation be- 
tween the FR and the emission line ratio. The black line in the 
plot corresponds to a least-squares-fit / = 51x 043 with the 
emission line ratio x. The [CII] emission is strongest for low 
to intermediate densities while HCO + is a typical density tracer. 
Both lines are typical PDR tracers, sensitive to significant FUV 
illumination. Consequently, a line ratio higher than 10 3 signals 
FR>ER. The corresponding points in the plot belong to models 
with n w 10 6 cm -3 and \ > 10 4 - The points with a FR<ER 
correspond to models with even higher densities and somewhat 
lower x- 

The panel on the right side shows the C + column density 
ratio plotted against T mb ( 13 CO(l - 0))/T mb ([CII]). All models 
with a line ratio higher than 2 contain fractionated C + . We saw 
in Fig. [10] that the FR(C + ) is strongest for low UV models with 
high densities. These models have weak C + emission, usually 
smaller than a few K km s" 1 , while these conditions favor 13 CO 
emission resulting in large line ratios. 

4. Summary 

We present an update of the isotope chemistry used in our PDR 
model code KOSMA-r. An automated routine was created to 
allow for the inclusion of isotope reactions into the chemical 
database files that are used in numerical PDR computations. This 
is combined with a proper rescaling of the new isotope reac- 
tions. We computed a large parameter grid of spherical PDR 
model clumps and investigated the effect of the isotope chem- 
istry, particularly that of the isotope exchange reaction ( IC U . on 
the chemical structure of the model clumps as well as on their 
emission characteristics. 

In the transition from ionized carbon to carbon monoxide the 
fractionation ratio of C + is always larger than the elemental ra- 
tio in the gas. Strong C + fractionation is possible in cool C + gas. 
However, this is only partly visible in the corresponding inten- 
sity ratios of the model clumps. Optical thickness and excitation 
effects produce intensity ratios between 40 and 400, strongly de- 
pendant on the model parameters. 

In the dense (n > 10 3 cm~ 3 ) gas, CO behaves differently 
and is never found with a fractionation ratio larger than the el- 
ement ratio with the exception of a very limited Ay ran ge un- 
der ve ry special parameter conditions. In the diffuse gas iLiszt] 
(2007) found qualitatively different behavior for n < 10 2 cm -3 



and M > 10 3 M Q . It turns out that isotope-selective photo- 
dissociation, the major process able to produce a FR>ER, is al- 
ways dominated by the chemistry in the denser PDR gas. This 
also affects the depth at which the transition from C + to CO oc- 
curs. The formation and destruction of 13 CO is much stronger 
controlled by reaction (IC U . A direct consequence is that in all 
models in our grid 13 CO is formed at smaller Ay than CO de- 
spite the weaker shielding capabilities of the rarer isotopologue. 

The fractionation of other species can be understood in terms 
of their formation history. If their major formation channel orig- 
inates from C + their FR is related to FR(C + ). This is the case for 
many light hydrides, especially CH, CHj, and CH3 . If the FR of 
the parental species is controlled by other reactions than ( IC Q 
the behavior might change. Atomic carbon is a mixed case with 
a regime at lower Ay where formation occurs mainly via recom- 
bination of C + and a regime where it is chemically derived from 
CO and HCO + . Consequently, the FR(C) exhibits a mixture of 
both cases. At particular depths, CH + might be formed from C. 
If that is the case the relation of its FR to FR(C + ) breaks down 
until formation via a C + route takes over again. 

Our computations have shown that CH and HCO + may be 
very sensitive tracers for chemical fractionation in PDRs. CH 
amplifies the FR of C + and HCO + amplifies the FR from CO in 
the CT. Both are abundant in the region where chemical fraction- 
ation plays a big role and suffer less from optical depth effects 
than their chemical progenitors. 

We demonstrated that fractionation of the local densities not 
necessarily transforms into a fractionated column density. C + is 
a prominent example. Fractionation of ionized carbon only takes 
part in cool C + gas which in most clouds only makes up a small 
fraction of the total gas column. Only low FUV models are able 
to produce larger columns of cool C + and will have a fraction- 
ated column density. CO behaves oppositely in that it requires 
warm CO in order to become fractionated. Again, this affects 
the column density only under certain conditions. 

Finally, we provide diagnostics for the fractionation status of 
C + , C and CO through suitable emission line ratios. We showed 
that a line ratio of T mb ( 13 CO(2 - l))/T mb ([CI] 610/im ) < 1 
signals a significant fractionation of the CO column density. The 
line ratio of T mb ([CI] 610Atm )/T mb (H 13 CO+(l-0)) has apower 
law dependence on the fractionation ratio of the column density 
of atomic carbon. The column density fractionation of ionized 
carbon is reflected in T mb ( 13 CO(l - 0))/T mb ([CII]) > 2. 
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Appendix A: Isotopization rules 

Usually, only the main isotopes are considered in astrochemical 
databases, despite the fact that many isotopologues have been 
detected in astronomical observations so far. We will describe 
here our efforts to include I3 C and I8 isotopes into the chem- 
ical database. However, in this paper we will only discuss the 
scientific implications of carbon fractionation not going into the 
details of the ls O chemistry. 

Usually, in a reaction formula like H + 2 — > OH + O it 
is not possible to identify which oxygen atom binds with hydro- 
gen. Both atoms in 2 are indistinguishable. This changes by 
including an isotope into a molecule. Now both atoms in 18 
are distinguishable and we get: 



H + 18 — > OH + 
H + 18 — ► 18 OH 



18 Q 

+ 



Unfortunately, this becomes more complicated if the isotope 
can be placed in more than one spot. The main isotope reaction 



C + HCO+ - 
splits into 5 isotopic reactions 

C + HC 18 0+ - 
C + H 13 CO+ - 

c + h 13 c 18q+ . 

13 C + HCO+ - 
13 C + HC 18 0+ - 



CO + CH 



C 18 + CH+ 

13 C0 + CH+ 
i3 C i8 + CH + 

CO + 13 CH+ 
C 18 + 13 CH+ 



taking into account that the C=0 binding is preserved (again 
we omit the isotopic superscript when denoting the main iso- 
tope). For more complex species the above scheme becomes 
much more complicated. Furthermore, the number of additional 
reactions is so large that it becomes difficult to perform the iso- 
topization by hand, especially if one plans to update the chem- 
ical network regularly and manual isotopization is quite error- 
prone. We developed a software routine to automatically imple- 
ment isotopic reactions into a given reaction se{3 The routine 
features are: 

- inclusion of a single 13 C and a single I8 isotope (multiple 
isotopizations are neglected in this study) 

- UDfA often does not give structural information, for in- 
stance C 2 H 3 does not distinguish between linear and circu- 
lar configurations (1— C 2 H 3 and c— C 2 H 3 ). In such cases we 
consider all carbon atoms (denoted by C n ) as indistinguish- 
ablfZli.e. 



However, if structure information is provided we account for 
each possible isotopologue individually: 



C 


+ H 2 CCC - 


-> C 4 H + H 




c + 


H 2 CC 13 C - 


-> C 3 13 CH-f 


H 


c + 


H 2 C 13 CC - 


-> C 3 13 CH-f 


H 


c + 


H 2 13 CCC - 


-> C 3 13 CH-f 


H 


13 C 


+ H 2 CCC - 


-> C 3 13 CH-f 


H 



molecular symmetries are preserved, i.e. NC CN = 
N 13 CCN, but HC 18 OOH ^ HC0 18 OH 
functional groups lik e CH n are preserved (see also 
IWoods & Willacvf 2009). e.g.: 



CH 2 CN H 
CH 9 13 CN H 



13 



but not: 



CHIN- 



CH, 13 CN + + e~ 



+ e~ 
+ e~ 
+ e~ 



CN + 
13 CN 
CN + 

CN- 



CH 9 



13 



CH 2 
CH 9 



13 



CH 9 



- when the above assumptions are in conflict to each other we 
assume minimal scrambling, i.e. we choose reactions such, 
that the fewest possible number of particles switch partners. 
For example: 



CH 3 OH -+ 


C 3 H+ - 


-> HC 3 0+ + CH 4 


CH 3 OH + C 2 


13 CH+ - 


->HC 2 13 CO+ + CH 4 


CH 3 18 OH4 


C 3 H+ - 


^HC 3 18 0+ + CH 4 


CH 3 18 OH + C 2 


13 CH+ - 


^HC 2 13 C 18 0+ + CH 4 


13 CH 3 OH + 


C 3 H+- 


-> HC 3 0+ + 13 CH 4 


13 CH 3 18 OH-F 


C 3 H+ - 


^HC 3 18 0+ + 13 CH 4 



but not 



13 



CH 3 OH 



C 3 H- 



HC 9 13 CO 



CH 4 



which would preserve the 3 C=0 bond but would require 5 
particles to switch partners. 

- we favor proton/H transfer over transfer of heavier atoms 

- we favor destruction of weaker bonds, e.g.: 



C+ + ocs 
c+ + o 13 cs 
C+ + 18 ocs 
c+ + 18 o 13 cs 
13 c+ + ocs 

13 C + + 18 0CS 



CS+ + 
CS+ + 
CS+ + 
CS+ + 

13 cs + 
13 cs+ 



CO 

13 co 
c 18 o 

13 C 18 

+ C0 

+ c 18 o 



C + C 3 H 2 
C + C 2 13 CH 2 



13 



C 4 H 



H 



C 3 13 CH 



C + C 3 H 2 — > C 3 CH 



H 
H. 



16 We realized the isotopization routine in Mathematica©by Wolfram 
Research. 

17 Notation remark: isotopization of C„ leads to C^/'C , e.g. 
C 2 > C 13 C. We always quote the full set of new reaction in the iso- 
topic network derived from one UDfA reaction 



but not: 



C+ + o 13 cs — ► 13 CS+ + CO 



In the above example the binding enthalpies in 0=C=S are 745 
kJ/mol for 0=C and 536kJ/mol C=S. We only allow reactions, 
that break the weaker bond. In other words, we preserve 0=C 
bounds above others, since it is the strongest double bound. 
Please note that this is a strong assumption that can be sacri- 
ficed if more detailed knowledge on the reaction kinetics is at 
hand. 
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Appendix B: Rescaling of reaction rates 

By introducing isotopologues into the chemistry we introduce 
many new reaction channels and we need to make sure that re- 
action rates are properly scaled. Unfortunately, not only the re- 
action rates for isotopologue reactions are unknown, we neither 
have information on branching ratios for reactions with several 
possible product channels. Thus, we assume equal probabilities 
for all branches, i.e. all isotopologue reactions possessing the 
same reactants but different products have their rate coefficient 
divided by the number of different product branches. For exam- 
ple: 



C + 9 



co + o 



4.7(-ll) 



with a being the fit coefficient from UDfA06 at T = 300 K 
and the number in parentheses indicates the decimal power. 
Introducing 13 C and ls O into this reactions opens up additional 
channels: 



C 
C 

13 

13 c 

13 C . 



18 
18 

c + o 2 

18 
18 



CO + 

c 18 o 
13 co 
13 co 

13q18 



o 



o 
o 

18q 

+ 



a = 2.35(-ll) 
a = 2.35(-ll) 
a = 4.7(-ll) 
a = 2.35(-ll) 
a = 2.35(-ll) 



Assuming equal probabilities for different reaction branches is 
a strong assumption and the reader should keep in mind, that 
many of the introduced reactions might have different reaction 
rate coefficients with potentially strong impact on the solution of 
the chemical network. The isotope exchange reactions discussed 
in the following are a notable exception from this. 

Appendix C: Influence of chemical data sets 

To illustrate how the choice of the chemical data set affects 
the outcome of our astrochemical calculations we calculate 
our reference model for three different chemical sets: UDfA06 
dWoodall et al.ll2007l) . O SU (version osu_01 2009), and KIDA 
(version kida.uva. 2011, Wakelametal. 2012). To allow a con- 
sistent computation it was necessary to slightly alter the original 
data sets: 

- The formation of H 2 on grain surfaces is calculated sep- 
arately in KOSMA-r and we removed the reaction H + 



H 



H 2 from OSU. 



- The same holds for the photo-dissociation of H 2 , we re- 
moved the reaction H 2 +h^ — > H+H. KOSMA-r explicitly 
calculates the H 2 formation rate from the population of the 
(vibrational) w-levels (all rotational levels of one vibrational 
state summed, ground state v = — 14, Lyman band 24 level 
+ Werner band 10 level) 

- The unshielded CO photo-dissociation rate coefficient dif- 
fers significantly among the three sets. We recomputed the 
unshielded photo-dissociation rate for a standard Draine 
FUV field using absorption cross sections available from 
http ://home . strw. leidenuniv.nl/^ewine/phot o/ and can con- 
firm the value of 2 x 10~ iU s -1 ( Ivan Dishoeck et al J 120061) 
from UDfA06. We replaced the corresponding a values 
in OSU and KIDA. F or a discussion on differing photo- 
rea ction rates see e.g. Ivan Hemert & van Dishoec k (2008) 
and lRollig et all (120121) . 

- We replaced reactions with very large negative rat e coef- 
ficient s 7 with a refitted expressions according to iRolligj 
(120111) . 



KIDA uses the Su-Chesnavich capture approach to com- 
pute rate coefficients for unmeasured reactions be tween ions 
and n eutral species with a dipole moment dWoon & Herbstl 
2009). Their formalism is incompatible with KOSMA-r and 
we replace the corresponding 1877 reaction rate coefficients 
with a new set of rate coefficients a, (3, and 7, suitable 
for the Arrhenius-Kooij formula a(T/300i^) /3 exp(-7T), 
which are fitted such that they approximate the original rates 
between 10 and 1000 K. 

For species where KOSMA-r does not distinguish between 
linear and cyclic isomeric forms we consider all species only 
in terms of their molecular formula, e.g. C 3 H instead of 
1 - C 3 H and c - C 3 H. 




Fig. C.l. Chemical structure of a model clump with the fol- 
lowing model parameters: no = 10 5 cm~ 3 , M = 100 Mq, 
X = 10. (UDfA06 (solid), OSU (dashed), KIDA (dotted). The 
gray shaded area shows the gas temperature spanned by the three 
models calculations. 




A v 

Fig. C.l. Same as Fig.|CJ]for HCO + , CH, and CH + . 

In Fig. IC.ll we show the relative abundances of C + , C, and 
CO for model calculations with n = 10 5 cm~ 3 , M = 1000 M©, 
and \ = 10 using different chemical data sets: UDfA06 (solid 
lines), OSU (dashed lines), KIDA (dotted lines) and the resulting 
range of gas temperatures. All three sets give the same chemical 
structure but with slight variations in the abundance profiles. The 
biggest difference is a higher CO abundance of OSU and KIDA 
at very low values of Ay (factor 2-3 at Ay = 0.01). This also 
leads to a carbon transition from C + to CO at slightly lower val- 
ues of Ay compared to UDfA. Overall, the chemical structure 
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of the main carbon species is comparable, particularly the effect 
on the total column density is small. The similar chemical struc- 
ture and the almost identical gas temperatures around the carbon 
transition will result in a consistent fractionation behavior of the 
three species across all three chemical sets. 

Fig. lC.2l shows how the different chemical sets influence the 
chemistry of HCO + , CH, and CH + . The differences are larger 
than those in Fig. IC.ll UDfA06 produces a significantly higher 
CH abundance in the outer layer of the cloud, at some positions 
by a factor of 10 higher than OSU. However, we discussed ear- 
lier that fractionation of CH is the result of C + fractionation. We 
do not expect any significant deviations in the fractionation ra- 
tio (FR) of C + ) for the three chemical sets, hence the effect on 
the FR(CH) should be weak. The same should be the case for 
CH + which shows even smaller differences for the three chem- 
ical sets. HCO + shows the same abundance in the cloud center, 
but some significant differences for Ay < 1. This could lead 
to different fractionation behavior at these parts of the cloud 
when using different chemical sets. It is unlikely that this leads to 
observable effects because of the 100 times lower HCO + abun- 
dance with respect to the center of the cloud and the respective 
weak influence on the total column density. 

Appendix D: Spherical model context 

D. 1. Cloud radius 

The model primarily computes radius dependent quantities. For 
a given density law, i.e., values a and f c = R C me/ Rtot, the total 
cloud radius i? t ot (in cm) can be calculated from uq (in cm -3 ) 
and M (in M©) using 



R tnt = 6.57 x 10 1 



(a - 3)M/° 



cm 



For a 



V "o (3/° - a/3) 
3/2 and f c = 0.2 this reduces to i?tot 



(D.l) 

5.3 x 



10 18 s/M/n cm. 



D.2. Column density 

Thegeneral expression for the maximum (radial) column den- 
sit)0is 



n Rtot (a f, 



l-a 



-1 



a-l 



that is iV max = A.lnR cm" 2 for a = 3/2 and f c = 0.2. 

However, they are not observables. Observations always 
yield an projected, beam-convolved figure. We describe them in 
terms of measurable column densities. 

In the framework of spherical model clouds, column densi- 
ties differ depending on where we look at. To get a position- 
independent measure for the column density of a given species 
i, we calculate the average column density for the whole clump 



<Ni > -- 



An 
^R 2 



rii(r)r 2 dr 



(D.3) 



When referring to model column densities we always mean a 
clump averaged column density according to Eq. ID. 31 This def- 
inition assumes that observations always cover whole clumps. 
This is equivalent to the traditional description of beam-filling 



JV max is the total gas column density from the cloud surface to the 
center along the cloud radius. 



factors for the observation of spatially unresolved clumps, i.e. we 
assume that we have (typically many) unresolved clumps within 
the telescope beam. In that sense, the fractionation of the column 
density is always the result of a convolution of the fractionation 
structure with the absolute abundance profile. 

D.3. CH column densities 

Modelling the formation of the light hydrides, such as CH + and 
CH, still poses a challenge to chemical models. In Fig. IC.3I we 
plot the mean CH column density versus the total mean column 
density of the respective clump. The coloured dots are observa- 
tions from absorption and emission line measurements. The red 
and black lines a re simple parametrized models to describe the 
column density (Mattila 19861: lOin et al.ll2010t) . We note, that 
the column density splits into two distinct regimes with signifi- 
cantly different behavior. Diffuse clouds with total column den- 
sities below 10 21 cm~ 2 show a steeper slope than denser clouds, 
where the CH column density appears to approach a limit of 
about 10 15 cm~ 2 . This behavior is approximately reproduced 
by the model calculations. The column densitie s in our model 
are co nsistent with values for TMC-1 given by ISuutarinen et ail 
(1201 lh as well as with diffuse cloud observations presented by 
iGerin et all d2010h which give JV(CH) w 1 - 26 x 10 13 cm~ 2 
for total H ? columns be tween 10 — 10 22 cm~ 2 . 

ISheffer et al.1 (1200 8) present column densities along 42 dif- 
fuse molecular Galactic sight lines. They find total columns be- 
tween 10 12 and 10 14 cm -2 and a v ery strong corre lation between 
column densities of CH and Ho. iMattilal ( 11986) confirms this 
trend for dark clouds. Magna ni et alJ ((2005) derived CH column 
densities from 3335 MHz observa tions in t he Gal actic plane and 
used the linear relation given bv iMatt ila (1986) to derive cor- 
responding H 2 column densities. They find 10 13 < iV(CH) < 

10 1 5 cm" 2 . 

IMagnani et all d2006h observed the CH 2 U m , J — 1/2, F = 
1 — 1 transition toward the Galactic center. They find iV(CH) « 
3 — 7 x 10 15 cm -2 . In addition to determining N(H 2 ) from its 
linear correlation to iV(CH) they also used a factor 1.8 x 10 20 
to derive jV(H 9 ) fro m integrated CO(1-0) line emission. 

lOin et all (HoiO) presented recent Herschel/HIFI observa- 
tions against Sgr B2(M) which revealed that the linear rela- 
tionship between CH and H 2 flattens at higher visual extinc- 
tions. They give a log-log slope of 0.38±0.07 for N(H 2 ) > 



(D.2) 10 21 cm 2 . In Fig. IC.3I we show our model results for 
iV(CH) versus N(R ) together with the observational data 
dRachford et alpOOl IMagnani et alj|200l 120061 ISheffer et all 
l2008blOin et alfcOlOLlGerin et al]l2010tlSuutarinen et alJl201 ll) . 

We cannot reprod uce the large column densities derived by 
Mag nani et al.l d2006l) . Their data is derived from emission line 
measurements, while all other studies used absorption lines, and 
consequently suffers from higher uncertainties due to its depen- 
dence on the assumed excitation temperature. To our knowl- 
edge 13 CH has not yet been detected. The lack of 13 CH obser- 
vations makes it difficult to assess the model results. The total 
CH column densities that we find in our model results stay be- 
low 10 15 cm -2 . Nevertheless, the general behavior is well re- 
produced. 



Appendix E: Fractionation plots of selected 
species 

Here we show fractionation plots of selected species over 
a large portion of our model parameter space. We left 




Fig. C.3. Mean CH col umn density v e rsus H 7 colum n density (color / shape coding same as F ig.|9]i. Observational results are plotte d 
as colo red points (red : ISheffer et al.1 d2008h: green: iRachford et al.1 (120021); olive: iMaenani et al.l (l2005h : black: lOin et all (120101) : 

blue lSuutarinen et all (1201 ll) ). The two lines are the CH-H 2 relations from 



brown : Magnani et al. (2006); orange Gerin et al.l d2Ql(J 
iMattilal (119861) (redl and loin et all d2oToh ("blackl. 



out models with M = 10 3 M Q because they show 
very similar results as models with M = 10 2 M . 
The fractionation plots are also available for download: 

|http : / /www . astro . uni-koeln . de/kosma-tau| 
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Fig. E.l. C + fractionation structure as function of relative clump radius r/R to t for different values of n and M. Each sector corre- 
sponds to a different \ value. The FR is color coded, ratios within ±10% of the ER are shown in green. Blue and violet denotes 
FR<ER, yellow and red denotes FR>ER. 




.E.2. Same as FiglETTIfor C. 
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Fig.E.3. Same as FiglRllfor CO. 




Fig.E.4. Same as FiglRllfor HCO + . 




.E.5. Same as FiglETTIfor CH. 




Fig.E.6. Same as FiglRllfor CH + . 



